
source("R/fct_Orthant.R")
source("R/Data_Orthant.R")
t<-proc.time()
reso<-order(A,a,m,M)
tOrd<-(proc.time()-t)

t<-proc.time()
resj<-jt(A,a,m,M,1)
j<-(proc.time()-t)


#t<-proc.time()
#resist<-ist(A,a,m,M)
#is<-(proc.time()-t)

t<-proc.time()
flt1<-flt(reso$S,reso$a,m,M,1)
fl<-(proc.time()-t)+tOrd
	
t<-proc.time()
b<-jt(reso$S,reso$a,m,M,1)
jj<-(proc.time()-t)+tOrd

t<-proc.time()
b<-jt(reso$S,reso$a,m,M,2)
j_RD<-(proc.time()-t)+tOrd

t<-proc.time()
res<-kt(A,a,m,M)
k<-(proc.time()-t)

t<-proc.time()
res<-kt(reso$S,reso$a,m,M)
kk<-(proc.time()-t)+tOrd


t<-proc.time()
#res<-it(A,a,m,M)
i<-proc.time()-t

t<-proc.time()
solm<-mt(reso$S,reso$a,m,M,2)
mi<-proc.time()-t+tOrd

solmst<-mst(reso$S,reso$a,m,M)
#res2<-gt(A,a,m)

cat("\n")
cat("\n")
t<-proc.time()
res3<-kt(A,a,m,M)
h<-(proc.time()-t)

Mh<-floor(j[[3]]/h[[3]]*M)
Mk<-floor(j[[3]]/k[[3]]*M)
Mjj<-floor(j[[3]]/jj[[3]]*M)
Mj_RD<-floor(j[[3]]/j_RD[[3]]*M)
Mmi<-floor(j[[3]]/mi[[3]]*M)
Mkk<-floor(j[[3]]/kk[[3]]*M)
Mfl<-floor(j[[3]]/fl[[3]]*M)
#Mist<-floor(j[[3]]/is[[3]]*M)
#Mj<-floor(f[[3]]/j[[3]]*M)

cat("\n")
cat("\n")
cat(resj)
cat("\n")
cat(res3)
#source("R/test_Orthan_S.R")
#solst<-jst(A,a,m,M/10,1,1,10)
#
Mm<-10
#XX<-rtmvnorm(10000,mean=m,sigma=A,lower=a,upper=rep(Inf,p))
i=1
b<-matrix(-Inf,nrow=Mm,ncol=7)
for(i in 1:Mm)
{
	b[i,1]<-jt(A,a,m,M,1)
#	b[i,2]<-(kt(A,a,m,Mk))
	b[i,3]<-jt(reso$S,reso$a,m,Mjj,1)
#	b[i,4]<-kt(reso$S,reso$a,m,Mkk)
	b[i,5]<-jt(reso$S,reso$a,m,Mj_RD,2)
	b[i,6]<-mt(reso$S,reso$a,m,Mmi,2)
	b[i,7]<-flt(reso$S,reso$a,m,Mjj,1)
#	b[i,7]<-(ist(A,a,m,Mist))
	
	#write.table(b,"res_Orthan.txt")
	cat("loop ")
	cat(i)
	cat("\n")
}
#.C("testbrt")
z<-factor(c(rep("Res",Mm),rep("NRes",Mm),rep("IRes",Mm),rep("INRes",Mm),rep("I-Res_OR",Mm),rep("BS",Mm), rep("FL",Mm)))
#z<-factor(c(rep("No-Res",Mm),rep("I-Res",Mm),rep("I-NoRes",Mm),rep("I-Res_RD",Mm)))
boxplot(as.vector(exp(b[,1:7]))~z)

pdf("~/Dropbox/ResultsGhk.pdf")
boxplot(as.vector(exp(b[,1:7]))~z)
dev.off()
dyn.unload("test/Orthant.so")
#dyn.unload("test/RLogit.so")
#YY<-rmvnorm(10000,mean=apply(XX,2,mean),sigma=var(XX))
#plot(YY)
#points(XX,col="red")
#dyn.unload("test/Rprobit.so")
#dyn.unload("test/IBIS.so")
